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We study the unpolarized fragmentation functions and parton distribution functions of the pion 
employing the nonlocal chiral quark model. This model manifests the nonlocal interactions between 
the quarks and pseudoscalar mesons in the light-cone coordinate. It turns out that the nonlocal in- 
teractions result in substantial differences in comparison to typical models with only local couplings. 
We also perform the high Q^-evolution for our results calculated at a relatively low renormalization 
scale ^ 0.36 GeV^ to be compared with the experimental data. Our results after evolution are 
in qualitatively good agreement with those data. 
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I. INTRODUCTION 

To apply perturbative Quantum Chromo dynamics (QCD) to study hadronic processes, such as the deep-inelastic 
electron-proton scattering, one needs the QCD factorization theorem to guarantee the cross section to be the convolu- 
tion of the two parts: the process-dependent perturbative QCD (pQCD) calculable short-distance parton cross section, 
and the universal long-distance functions which can be extracted from experiments. Since the strong interaction in 
long-distance is essentially nonperturbative, those long-distance functions are incalculable by pQCD. 

To date, there has been increasing interest in those long-distance functions such as the fragmentation functions 
and parton distribution functions. The fragmentation function characterizes the probability for a hadron fragmented 
from a quark with the momentum fraction z. It plays a crucial role in analyzing the semi-inclusive processes in 
the electron-positron scattering, deep-inelastic proton-proton scattering, and so on [THll]. The parton distribution 
function provides information for the distribution of the momentum fraction x carried by a parton inside a hadron. 
It is a necessary ingredient in studying hard scattering processes such as the deep-inelastic electron-proton scattering. 
It is worth mentioning that a parton distribution function can be extracted from the corresponding fragmentation 
function via the Drell-Levy-Yan (DLY) relation due to the analytical continuation [12 . This analytic continuation is 
possible only if the parton distribution function and the fragmentation function are assumed to be described by the 
same function defined in the different regions [12 . 

Because of these salient physical meanings implied in those functions, they have been intensively studied so far ^l3l - 
[38] . These functions have been extracted from the available empirical data by global analyses or from the parametriza- 
tions which satisfy certain constraints p^SHTS] . Many theoretical investigations on the fragmentation functions have 
been carried out. For example, a rigorous QCD proof has done for the momentum sum rules of the fragmentation 
functions [21]. The Dyson-Schwinger (DS) method was applied to compute the valance-quark distribution function, 
resulting in relatively good agreement with the available experimental data for the pion and kaon [27^. Similarly, 
the pion parton distribution was studied with a nonlocal Lagrangian using the DS method [28]. In Ref. [29], the 
valance-quark distribution was determined by analyzing the Drell-Yan process including the soft-gluon resummation. 
A review for the experimental and theoretical status for valance-quark distribution for the nucleon and pion is given 
in Ref. [30 . The authors of Refs. j3TH34] have also made use of the instanton- motivated approaches to compute 
the quark distribution functions as well. The collins fragmentation functions which play an important role in the 
transverse spin physics have also been studied [35H37] . Monte-Carlo simulations with supersymmetric (SUSY) QCD 
were carried out to obtain the fragmentation function up to a very high energy in the center-of-mass frame y/s [38] . 

Among those approaches, we are in particular interested in the ones based on the chiral dynamics. Because the 
strong interaction in long-distance is dominated by the chiral physics, therefore the nonperturbative objects such as 
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fragmentation functions also should be subjected to the chiral physics. Along this line, the fragmentation functions 
were first studied by Georgi-Manohar model in Refs. [Hill], and later by chiral-quark- meson coupling models in the 
pseudoscalar (PS) and pseudovector (PV) schemes [20l [35l [36] . Their results indicate considerable differences between 
these two schemes. In Refs. [22H26] , the authors studied the fragmentation functions for various hadrons in terms 
of the Nambu-Jona-Lasinio (NJL) model. In their works, quark-jets and resonances were taken into account and 
consequently the momentum sum rules is satisfied. The quark-jet contributions turned out to be crucial to reproduce 
the various fragmentation functions for the small z region. Besides, the quark-jet is also essential to generate the 
unfavored fragmentation functions [23, 25, 26 . 

In the present work, we concentrate on the fragmentation functions and quark distribution functions for the SU(2) 
light flavor for the positively charged pion, i.e. 7r+ = ud. We only evaluate the 7r+ case because the other isospin 
channels can be easily obtain by multiplying the corresponding isospin factors in our model [25]. For this purpose, 
we employ the nonlocal chiral quark model (NLChQM). This model is based on the dilute instanton-liquid model 
(LIM), which is properly defined in Euclidean space [4Q.-;44 . It is worth mentioning that, NLChQM have successfully 
produced results for various nonperturbative QCD quantities which are in agreement with experiments as well as 
lattice QCD (LQCD) simulations. [45H5Q) . 

Our calculations are performed in the light-cone coordinate at a relatively low renormalization scale \i ~ 0.6 GeV, 
corresponding to the two phenomenological instanton parameters, i.e. average size of the (anti)instanton p ^ 1/3 fm 
and average inter-instanton distance ^ ^ 1 fm for the dilute (ant i) instanton ensemble. These instanton parameters 
have been determined phenomenologically [40 or generated by LQCD simulations [51 . Although the fragmentation 
functions and the parton distribution functions are defined properly only in Minkowski space and there has been no 
rigorous proof, we still assume that there is an appropriate analytic continuation between the instanton physics and 
the NLChQM defined in the light-cone coordinate. We notice that there have been successful applications based on 
this assumption to various physical quantities (311434] . 

Distinctive features of the present approach, in comparison to other local effective chiral models such as the con- 
ventional NJL model and usual quark-meson coupling models, are as follows: 1) The couplings between quarks and 
pseudoscalar (PS) mesons are nonlocal. Those nonlocal interactions are led by the intricate quark- (anti) instanton 
interactions inside the instanton ensemble, flipping quark helicities [4QH44] . 2) As a result, the quark-PS meson 
coupling, which is represented by the momentum-dependent constituent- quark mass, plays the role of the natural 
UV regulator. Hence, one does not need any artificial regulators. Interestingly enough, the computed momentum- 
dependence of the mass turns out to be very comparable with the LQCD simulations [^|^. 3) Moreover, all the 
relevant physical quantities such as the low-energy constants (LECs) are determined self-consistently only with the 
instanton parameters [53-55 . There is no more adjustable parameters in principle within the framework as long as 
we are considering the light-flavor SU(2) sector. We emphasize that, among the above features, the nonlocality of the 
interaction is the most important ingredient which has never appeared in other chiral models. 

As for the numerical calculations, we compute the elementary fragmentation function d^{z\ standing for the 
fragmentation process u nq^ where tt is assigned for 7r+ throughout this work, i.e. the quark u quark is fragmented 
into the pion and an energetic quark q\ Since an energetic quark will be fragmented into other hadrons until it 
possesses insufficient energy to be fragmented into hadrons, this elementary fragmentation function does not satisfy 
the momentum sum rule 231125]. It is likely that if the quark jets and resonance contributions are taken into account, 
the sum rule will be satisfied [23] [25]. However we will not consider those effects for brevity in the present work. Our 
results of the fragmentation functions turn out to be substantially larger in comparison to other models. Actually the 
first moment of the elementary fragmentation function can be express as [23] 




Here Zq is the residue of the quark propagator in the presence of the pion cloud and can be interpreted as the 
probability of finding a bare constituent quark without the pion cloud. In other words, those elementary fragmentation 
functions are normalized to the number of pions per quark ^23^ . The value of Zq in our model is much smaller compared 
with other models with local couplings. Physically we can understand the difference between our model and other 
models is due to that the nonlocal couplings between quarks and PS mesons in our model somehow mimic a part 
of the meson-cloud effects. Therefore our elementary fragmentation functions have already contained more pion- 
cloud contributions than other models. To investigate the z-dependence of the fragmentation functions from different 
models, we renormalize these elementary fragmentation functions as follows 
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FIG. 1: Schematic figures for tfie quark-distribution (left) and fragmentation (right) functions, in which the soUd and dash 
hues denote the quark and pseudoscalar meson, respectively. 

We find that our fragmentation functions multiplied by i.e. zD^^ is more symmetric with respect to z compared 
with the ones from other models. 

The calculated fragmentation functions are evolved to high by the DGLAP evolution using the code QCD- 
NUM |60l [61] to compare with the empirical data. By doing so, we see qualitative agreement between our result 
and the empirical data, although some underestimates are shown in the small-z region. This defect is expected to be 
improved by the inclusion of the quark-jet contributions [23[ 125] . 

The parton distribution function f^{x) is extracted from the fragmentation function by the DLY relation. As 
done in Ref. [23], we provide the numerical results for the minus- type and plus- type quark distribution functions 
being multiplied by x, i.e. x{f^ — f^) and x{f^ + /^), respectively. Note that the minus-type one is nothing but 
the valance quark-distribution function. Again, we perform the high evolution in order to compare our results 
with the empirical ones, resulting in qualitative agreement except the sizable underestimate observed again in the 
small-x region after the evolution up to = 4GeV^ for the minus-type one. Again, this discrepancy implies more 
realistic treatments beyond the elementary process, such as the quark-jet contributions as mentioned previously [23 , is 
required. However, implanting the quark-jet contributions into the present nonlocal-interaction models is challenging, 
since the momentum-dependent quark-PS meson coupling can provide nontrivial effects to the DLY relation and the 
quark-jet contributions here. So we leave it for our future study. It also turns out that the high evolution up to 
= 27 GeV^ reproduces the experimental data qualitatively well. 

The present report is structured as follows: In Section II, we briefly introduce the model we use and sketch our 
computation procedure. In Section III we present our numerical results and related discussions. The final Section is 
devoted for the summary, conclusion, and future perspectives. 



II. NONLOCAL PION-QUARK COUPLING 

The unpolarized fragmentation function indicates the process that an off-shell quark (q) is fragmented into an 
unobserved set of particles {X) and on-shell hadron (/i), i.e. q ^ hX. A schematic figure for this function is given in 
the right of Figure [ij Here, the fragmentation function is defined in the light-cone coordinate, assuming the light-cone 
gauge, as follows p5l l36]: 

D^iz.e^,^) = j^l d/c+Tr [A(/c,p,/i)7-] U-=p-. (3) 

where /c, p, and z stand for the four-momenta of the initial quark and fragmented hadron, and the longitudinal 
momentum fraction possessed by the hadron, respectively. k± denotes (ko ± ks)/V^ in the light-cone coordinate. 
All the calculations are carried out in the frame where k± = 0. Here the z-axis is chosen to be the direction of k. 
On the other hand, kr = k — (^^^^ Pi defined as the the transverse momentum of the initial quark with respect to 

the direction of the momentum of the produced hadron, is nonzero, /a denotes the renormalization scale at which 
the fragmentation process computed. Note that we consider this renormalization scale is almost the same with the 
momentum-transfer scale, i.e. /i^ ~ for simplicity, unless otherwise stated. The correlation A(/c,p, /i) reads 
generically: 

A(fc,p,M) = Y.J 0^e+^''<{om)\h,X){h,X\m\O)- (4) 

Here i/j represents the quark field, whereas ^ the spatial interval on the light cone. Furthermore one can integrate 
over /cx. 
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The factor of is due to the fact that the integration is over — V — (^p) transverse momentum of the 

produced hadron with respect to the quark direction, and there is a relation between and /ct- P± = —zkr The 
integrated fragmentation function satisfies the momentum sum rule: 



/ y2zD^{z,fi)dz = i, 



(6) 



where h indicates for all the possible hadrons fragmented. Eq. (|6| means that all of the momentum of the initial quark q 
is transferred into the momenta of the fragmented hadrons. From the Drell-Levi-Yan (DLY) relation p!2 ] [22l l23l [25| [36] , 
Dq can be related to the parton distribution function provided that there is a proper analytic continuation. The 
relation is as follows, 



(7) 



where x denotes the momentum fraction possessed by a parton inside the hadron. A schematic figure for the quark- 
distribution function is depicted in the left of Figure [l] 

In this article, we use NLChQM to investigate these nonperturbative objects, i.e. fragmentation and parton 
distribution functions. This model is motivated from the dilute instanton liquid model [4Qti44] . We note that, to 
date, various nonperturbative QCD properties have been well studied in terms of the instanton vacuum configuration 
and the results are comparable with experiments and LQCD simulations [45-48 . In that model, nonperturbative 
QCD effects are deciphered by the nontrivial quark-instanton interactions in the dilute instanton ensemble. However 
this model by nature is defined in Euclidean space because the (anti)instantons are well defined there by signaling 
the tunneling between the infinitely degenerate QCD vacua. Although there have been no rigorous derivation on the 
analytic continuation from the instanton physics to those in Minkowski one, there are still several challenging studies 
which try to apply the idea of the instanton physics to the physical quantities defined properly only in Minkowski 
space, such as the light-cone wave function j3TH34] . Following those studies we adopt the effective chiral action 
(EChA) from NLChQM in Minkowski space as follows: 



<Seff[m/,(/)] = -Spin 



(8) 



where Sp and rhf denote the functional trace Tr J d^x{x\ • • •\x) over all the relevant spin spaces and SU(2) current- 
quark mass matrix diag(m^i, m^i), respectively. Throughout the present work, we will take rriu = rrid = 5 MeV, 
considering the isospin space. Note that, in deriving EChA in Eq. Q, we simply replace the Euclidean metric for the 
(anti) instanton effective chiral action into that for Minkowski space. The momentum-dependent effective quark mass 
generated from the interactions between the quarks and nonperturbative QCD vacuum, can be written in a simple 
n-pole type form factor as follows 



M{d^) = Mo 



(9) 



where n indicates a positive integer number. We will take n = 2 as in the instanton model [33l [34]. A stands for 
the model renormalization scale. It is related to the average (anti) instanton size p. The nonlinear PS-meson field, 
i.e. U^^ takes a simple form [44] with the normalization following Ref. [35] to be consistent with the definition of the 
fragmentation function in Eq. ([3|: 



/7^^((^) =exp 



z75(t • 0) 
2F^ 



1 



2F^ 



(10) 



where stands for the weak-decay constant for the PS meson 0. For instance, is chosen to be about 93 MeV in 
this normalization. We note that, however, the value of F^^ can be determined rather phenomenologically to reproduce 
relevant physical quantities and conditions, even in NLChQM, and will discuss this in detail in Section III. We also 
write explicitly the PS-meson fields: 



T ■ (f) ■■ 



( "° 



V2' 



(11) 



By expanding the nonlinear PS-meson field up to (9(^1) from EChA in Eq. ([8|, one can derive the following effective 
interaction Lagrangian density in the coordinate space for the nonlocal quark-quark-PS meson vertex: 



-4>qq 



(12) 
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Then, the correlation in Eq. Q can be evaluated using this effective Lagrangian for the elementary fragmentation 
function q{k) (j){p) + q'{r)^ where r = k — p: 



1 MkMr 
1 MkMr 




75 {f + Mr) 75 

75 {f + M//) 75 



(2^)4 2F| 

where the momentum-dependent effective quark mass is written via Eq. (joj): 

2A2 




Me = Mo 



n 2 



2A2 - £2 + 



(13) 



(14) 



Note that the ja dependence of A(fe, p, /i) is due to that fact we need regularize Mi by introducing the model renor- 
malization scalar A in Eq. ([9|. In Eq. (13) we adopt the following notations: M^ = mj + Mi and Mf = rrif -\- Mq. 
Note that / and f indicate the flavors for the initial (q) and final (q^) quarks, respectively. Eq. (13) stands for the 
7r+ elementary fragmentation, i.e. u 7r+ d. As for the fragmentation of a neutral pion, one needs multiply the 
isospin factor 1/2 to Eq. (13) at the elementary-fragmentation level. From the first line to the second line in Eq. (13), 
we approximate some of the momentum- dependent effective quark masses into a constant constituent-quark mass 
Mi Mq, since the momentum dependencies from those mass terms play only minor roles numerically. Moreover, 
this approximation simplifies the numerical calculations to a great extent. Performing the trace over the Lorentz 
indices for the fragmentation function in Eq. ([s]): 



Tr [(^ + Mf) 76 (1^ - + M/') 76 + Mf) 7^] = -4 [p^{k'' - M]) + k^{M] - 2MfMf, + fc^ - 2fc • p)] , 
we reach a concise expression for the elementary fragmentation function u — > tt+cZ from NLChQM: 



MkMr 



87r3z(l-z) 2F2 



2(fc2 - Mf) + (fc2 + 



Mj 



2MfMf, -2k-p) 



(fc2 - M2)2 



(15) 



(16) 



where the relevant four momenta squared in the present theoretical calculations are defined in the light-cone coordinate 
as follows: 



e = k^k_ 



= {k-pf 



(kj^p- -^k_p^), 



where the transverse PS-meson momentum on the light cone, i.e. Pt = 0. Note that the value of Mq can be fixed 
self-consistent ly within the instanton model j4QH48] with the phenomenological (anti)instanton parameters p ~ 1/3 
fm and ^ ^ 1 fm, resulting in Mq ~ 350 MeV by the following self-consistent equation, defined in Euclidean {E) 
space: 



1 



d^q M{q^) 



M{q^) = Mo 



1 2 



p'^q'^ 



(17) 



The pion mass is chosen to be = 140 MeV throughout the present work. Collecting all the ingredients above, one 
is led to a final expression for the elementary fragmentation function: 



dl {z,kl,p) 



1 MkMr 



[z^kl^[{z-l)Mf^Mf>f] 



Stt^ 2F2 [z^kj.^z{z-l)M]^zM],^{l-z)ml]^'' 



(18) 



where M^ and Mr are the momentum-dependent quark mass manifesting the nonlocal quark-PS meson interactions, 
read: 



Mk = 



Mo[2A^z{l 



.)]2 



[^2fe2 + z{z - 1)(2A2 - ,52) + zMj, + {1- z)ml\^ ' 



Mr = 



Mo(2A2)2 

(2A2-M2,)2- 



(19) 



As for Mk in Eq. ( |19[ ), we have introduced a free and finite- valued parameter S in the denominator to avoid the 
unphysical singularities, which appear in the vicinity of (2:, /ct) = 0, due the present parametrization of the effective 
quark mass as in Eq. We determine 5 to satisfy a natural conditions for all the possible values of (2;, kr)' 



\z'^k^ + z{z ■ 



l){nh^ - J^) + zM]^ + (1 - z)ml\ > 0. 



(20) 
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It is easy to see that one of the trivial solutions for Eq. ( 20 ) can be 5'^ = n — Mj , considering the quark propagator 
in Eq. (18). Although there can be other solutions, we will use this value for 6 for the fragmentation function as a 
trial hereafter. Physically, those unphysical singularities can be understood as that a hypothetic particle, whose mass 
corresponds to a/2A, becomes its on-mass shell for a certain combination of (z, kr)- Thus, we employed S to exclude 
that unphysical situation. We also verified that the change of 6 does not make considerable effects on the numerical 
results, as far as Eq. (I20|) is fulfilled. 



If we replace all the momentum-dependent masses into a constant one M in Eqs. (18) and (19) in the chiral limit, 
and change the quark-PS meson coupling into a constant one appropriately in Eq. (18), we obtain the expression for 
the elementary fragmentation function from the quark-meson coupling model [36] : 



12- 



(21) 



The above equation is also equivalent with that from the N JL model calculations in principle [22l [23l [25] . Note that 
the on-sheh value of M^/{2F^) becomes about 7, which is quite similar to g^^^ ~ 9, used in Refs. [22l |23l [25l [36] . At 
the renormalization scale in our model, the elementary fragmentation function is assumed to be able to be evaluated 
further by integrating Eq. (18) over kr as: 



dl (2;,/i) 



2'Kz'^ I {z ^ ^ ji) Ut dkT ^ 

Jo 



(22) 



where the factor z'^ in the right-hand-side is again comes from the integration over p± = —kr/z. Actually the 
connection between and d^ {z^k'^^ji) would be far more complicated in principle [39] . We will present the 

numerical results for Eq. ( [22] ) in the next Section. 

Using the DLY relation in Eq. ([7|, one can derive the quark-distribution function for 7r+, i.e. /J^ as follows: 



3 MkMr 



fcf, + [{x - l)Mf - xMf.f 



47r3 2F2 [k'^ + {l-x)M'j + xM'j,+x{x-l)mlY' 
where we have defined the effective quark masses for the quark distribution function by 

4MoA4(1 -xf , ^ 4MoA4 



[fcf, + (1 - x){2h?) + xM'j, + x{x - l)m: 



212 ' 



Mr = 



(2A2 -M2^)2- 



(23) 



(24) 



Note that A1/c in Eq. ( |24[ ) does not suffer from the unphysical singularities unlike that in Eq. (19), according to 
the different kinematic situations between those functions. Again, if we take the replacements as Mjj/ M and 
A4kM.r/{2F^) g'^^^ as done for the fragmentation function, one is led to 



M2 



47^3 [fe2 + M2 + x{x - l)m2]2 ' 



(25) 



which is equivalent to Eq. (14) in Ref. [23 in the chiral limit. Similarly, the integration over can be performed as 
follows, resulting in a function of x at a certain renormalization scale /a ^ A: 



27r / {x,k^,lj)kTdkT. 
Jo 



(26) 



We will use the notation tt instead of 7r+ from now on for convenience. Moreover, the minus-type quark-distribution 
function, which corresponds to the valance-quark distribution, satisfies the following normalization condition: 



dx 



(27) 



where Uq in the right-hand-side indicates the valance-quark number and becomes unity. 



III. NUMERICAL RESULTS 



Here we present our numerical results for the fragmentation functions and the parton distribution functions of tt with 
relevant discussions. We will consider only the favored fragmentation process u nd. As in the instanton model the 
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value of A is proportional to 1/p giving the renormalization scale A ^ 600 MeV as mentioned in the previous section. 



As understood by Eq. (14), the effective quark mass plays the role of the UV regulator by construction, resulting in 
that relevant integrals in Eqs. ( [22| ) and ( [26| ) are not divergent. This feature is quite different from Refs. [22) l23l l25l [36] , 
in which various cutoff scheme were employed. Moreover, we assume that there is no perturbative gluon contributions 
for those function in hand at the initial scale A. 

First, we present the numerical results for the elementary fragmentation functions, multiplied by z, zd^ as functions 
of z in the left panel of Figure [2] for the present (solid), NJL (dot), PS (dash), and PV (dot-dash) results. As for 
the NJL result, the renormalization scale is chosen to be 0.18 GeV^ as usual [22j |23j [25l whereas it becomes 
1 GeV^ for the PS and PV schemes of the quark- meson coupling model [36]. It is worth mentioning that the initial 
values for each curve are chosen to be the same with each model renormalization scales. For instance, we set 
Q2 = (600 MeV) 2 for the present case. Although, in order for an appropriate comparison between the model results, 
one may need proper evolutions to a single value, we do not perform those evolutions here. Hence, the curves 
in Figure [2] should be taken into account as a very qualitative comparison at a generic hadron scale (0.4 ^ 1.0) GeV. 
From the left panel of Figure |2| we observe that our curves are in general much larger than the other model results. 
Moreover, while the peaks of the curves locate at z ~ 0.5 for the present theoretical framework, the maximum values 
of NJL and PS scheme are around z ~ 0.7. On the other hand, the maximum value occurs in PV scheme at the 
location where is close to ours. The behaviors of the NJL and PS-scheme curves are similar because they are both 
local models without derivatives. 

Note that these elementary fragmentation functions computed here do not satisfy the momentum sum rule Eq. (|6| . 
Since we have considered the elementary process q tt q\ however, as discussed in Refs. [23l[25], the pion-cloud 
effects will enhance the elementary fragmentation functions, and those fragmentation functions with the pion-cloud 
effects were identified as the renormalized fragmentation function. The pion-cloud effects can be estimated by the 
value of the first moment of the fragmentation functions: 

1-Zq = Y.I dzdl{z)=K/u- (28) 

where the summation runs over the all isospin states of the PS meson. For instance, in Refs. [23l[25], the value of JV^/u 
was estimated to be (0.1 ~ 0.2). As understood, this smallness of JV^/u results in the breakdown of the momentum 
sum rule. The renormalized fragmentation function then can be defined as [23^ : 

DZiz) = f^. (29) 

Here, we note that the summation runs over all the possible fragmented hadron states from the u quark in principle. 
However, since we are interested only in the pion case, we restrict ourselves to the SU(2) light-flavor mesons, i.e. no 
kaons. From the numerical computations, we obtain 

AP/^^^""* = 0.3685, AC/i^ = 0.0257, A/'J/^ = 0.0843, A/'J/^ = 0.1684. (30) 



From the values in Eq. (30), we find that present result gives much larger strength compared wit those in the models 
with local interactions. This observation tells us that the probability for the initial quark to be fragmented into the 
hadrons, especially pion, turns out to be much higher due to the nonlocal interactions. Physically, this enhancement 
of Afj^/q can be understood in the following way: The nonlocal interactions are originally generated from the intricate 
quark- (anti)instanton interactions in the instanton- vacuum picture [40-44 . This momentum-dependent quark mass 
can be understood as a dressed-quark mass as in usual Dyson-Schwinger methods [27] where the quarks are dressed by 
the pion cloud. In this sense, the nonlocal interaction represented by the momentum-dependent effective quark mass 
practically contains considerable contributions from the pion cloud. It leads to the higher probabilities for the initial 
quark to be fragmented into the pion, comparing to usual local models with a constant quark mass and coupling. The 
numerical results for the renormalized fragmentation function multiplied by z^ zD^{z) is given in the right panel of 
Figure [2] Interestingly, the curves from the present nonlocal model calculations look like an average of the PS- (or 
NJL) and PV-scheme results. This can be understood qualitatively by expanding the momentum-dependent quark 



masses of VM^M^ Mj, in Eq. (14): 



Mk^Mo-Mo^^--- . (31) 



The first term of Eq. (31) relates to the PS scheme (or NJL) without derivatives, while the second term to the PV 
one with derivatives, qualitatively. In this way, the present result can be seen as sort of a mixture of the two coupling 
models. 
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FIG. 2: (Color online) In the left panel, we show zdZ(z) at low renormalization scale for the present result (Q^ = 0.36 GeV^, 
solid), NJL 23, 25^ (Q^ = 0.18 GeV^ dot), PS scheme [35l|36] (Q^ = 1 GeV^ dash), and PV scheme [35l|3B] (Q^ = 1 GeV^ 
dot-dash). zDZ{z) using Eq. (29) are also shown in the right panel in the same manner with the left panel. 
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FIG. 3: (Color online) xf^{x) from the present model (solid), and NJL 



25] (dot) for ^ 93 MeV (left) and 85 MeV 



(right) at low values: = 0.36 GeV^ for the present models and = 0.18 GeV^ for the NJL. 



The numerical results for xf^ are given in Figure [s] for the present result (solid) and NJL (dot) cases as functions 
of the momentum fraction x. All the calculations are carried out at the same renormalization scale 600 MeV as for 



the fragmentation functions. If we integrate the quark-distribution function over x as in Eq. (27) with the empirical 
pion-decay constant F^^ ^ 93 MeV for Eq. (23), shown in the left panel of Figure |3j we have Uq = 0.88. These 
values deviate the normalization condition oiEq. (27) by (10 ~ 20)%. The reason for this discrepancy can be 



understood as follows: Within NLChQM (and also the instanton model), the value of F^ turns out to be about 
20% smaller than its empirical value, if the mass- derivative terms, i.e. dM}^/d\k\^ are ignored [4QH48] . Note that 
these terms are originated by considering the (axial) vector current conservation for the relevant hadronic matrix 
elements [411 |57l |58] . Practically, they are called as nonlocal contributions. In deriving in Section II, however, 
we did not include them and derived by employing the DLY relation. Consequently our values of violate the 
normalization condition Eq. (27). Thus, if we choose a smaller value for F^^, considering the absence of the nonlocal 
contributions, the normalization condition can be satisfied with F^^ ^ 85 MeV. In other words, the effects from the 
nonlocal contributions are compensated by the smaller F^ value. The inclusion of the nonlocal contributions is under 
progress [59 . The numerical results which satisfy Eq. (27) are given in the right panel of Figure [sj As shown there, 
the overall strength of the curves are enhanced by about 10%, whereas the curves shapes remain almost the same. 
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From a phenomenological point of view, we will use the numerical results for satisfying the normalization, given in 
the right panel of Figure [3j hereafter. It also turns out that the curve from the present calculation is very symmetric 
and peak at x ~ 0.5. The NJL result is tilted to the the region near x = 1. Thus, we conclude that the nonlocal 
interactions between the quark and PS meson have large impact on the shape of the curves of xf^ in comparison to 
those from the local-interaction model calculations, as already observed in the case of the fragmentation functions in 
Figure [2] 

According to Ref. [18 , xf^ can be parameterized by the following forms: 



x)^^ : quark contribution, 
: antiquark contribution. 



(32) 



where we have assumed that SU(2) flavor symmetry for the valance and sea quarks, respectively. Using Eq. (32), 
xf^{x)m Figure [3] can be parameterized at = 0.36 GeV^ and the fitted values for A^, 0^7^, and /^t^ are listed in 
Table llj satisfying the normalization condition in Eq. (27). These parameterized distribution functions will be used for 
the inputs for the high-Q^ evolution. To compare with the empirical data obtained in Ref. [14 , we need to perform the 
high-Q^ evolution of the present results. For this purpose, we will make use of the numerical DGLAP evolution code 
QCDNUM17 hy Botje [60l[6l]. The quark distribution functions computed and parameterized previously becomes 
the inputs for the evolution. For the QCDNUM evolution, one needs three valance quark x{q — q) and anti-quark xq 
distributions. Note that we assume that the gluon functions is zero at the initial scale. 



= 0.36 GeV^ 






P. 


Present 


3.3818 


1.2163 


1.4280 


NJL 


5.2497 


1.9131 


0.9208 



TABLE L Fitted values for the quark distribution function, AZ^ q^tt, and in Eq. (32) at Q = 0.36 GeV 



First, in Figure [ij we show the numerical results for the renormalized zD^ evolved from — 0.36 GeV^ (solid) to 
= 1 GeV^ at LO (dot) evolution. The empirical fragmentation function is parametrized as follows [M]: 



M7. 



(33) 



where B(x, y) in the denominator stands for the beta function with arguments x and y. The LO and NLO values for 
a and ^5, obtained at = 1 GeV are given in Table [n] following Ref. [M]. 





M: (LO) 


aZ (LO) 


K (LO) 


(NLO) 


aZ (NLO) 


(NLO) 


Dl(z) 


0.546 ± 0.085 


-1.100 ±0.183 


1.282 ±0.140 


0.401 ± 0.052 


-0.963 ±0.177 


1.370 ±0.144 



TABLE IL Parameters for the empirical pion fragmentation functions in Eq. (33), evaluated at = 1 GeV^ for LO and NLO 
evolutions, in Ref. 14 . 

The empirical results for zD^^ at = 1 GeV^ at LO and NLO are drawn in Figure [4] in the dash and long-dash 
lines, respectively with the inputs including Eq. (33) and the values listed in Table (llj The errors for the empirical 
fragmentation functions are also depicted. As for the present result (left), the LO and NLO empirical curves are 
qualitatively matching with the theoretical calculation (dot) for the region z > 0.5 whereas they deviate considerably 
in the region z < 0.5. As discussed in Refs. ^31 [25], the fragmentation functions will be enhanced by taking into 
account the quark-jet contributions in the region z < 0.5. Hence, although we have not taken those contributions 
into account here, the discrepancy between the theoretical and empirical curves might be cured by including the jet 
contributions which are now under progress [59] . 

We demonstrate our results evolved from = 0.36 GeV^ (solid) to 4 GeV^ at LO (dot) and NLO (dash) in Figu re 5 
for the minus- type (left column) and plus- type (right column). The empirical curves are computed using Eq. ( |32^ 
and input values fitted from the NAIO data [18] in Table IIL Note that we did not present the errors for the empirical 
data in Figure [sj since the errors given in Ref. [TF provide considerably small effects. It is clearly shown that the 
present results are in agreement with the empirical one for the minus-type distribution, which is nothing but the 
valance quark distribution function. We also observe that the present results for the plus-type distribution functions 
underestimate in the region x < 0.2 . 

After the evolution, we plot the valance quark distribution function at = 27GeV^ in Figure [o] The data points 
are taken from Ref. [16 . The present result reproduces the data qualitatively well as expected from the left panel of 
Figure [5] 



10 




FIG. 4: (Color online) Renormalized elementary fragmentation functions, multiplied by 2;, zD^^z) at — 0.36 GeV^ (solid) 
and they are evolved to — 1 GeV^ at LO (dot). The empirical data are taken from Ref. [14 for the LO (dash) and NLO 
(long-dash) analyses at = 1 GeV^. 







At 




0.64 ±0.03 


1.08 ±0.22 


0.9 ±0.3 


5.0 



TABLE III: Input values for the empirical quark distribution function for Eq. (32), given in Ref. "18". 



IV. SUMMARY AND CONCLUSION 



We have studied the fragmentation functions and the parton distribution functions for the positively charged pion 
using NLChQM motivated by the instanton vacuum configuration. We have computed them at the low-renormalization 
scale ^ ^ 0.36 GeV^, then evolved them up to certain higher values by the DGLAP evolution. The numerical 
results were compared with other models with only local interactions. Below, we list our important observations in 
this work: 

• The fragmentation and parton distribution functions calculated by NLChQM manifest the effects from the 
nonlocal-interactions between the quarks and the PS mesons. Naturally, if we reduce the nonlocal-interaction 
into local ones, we can obtain the well-known results such as those from the NJL model. The nonlocal interactions 
make significant differences in comparison to those from the models with local interactions. 

• By comparing the first moment of the elementary fragmentation function of all isospin states, Mj^/u from all 
models, the pion-cloud effects turns out to be much more pronounced in NLChQM since the intricate nonlocal 
quark-PS meson interaction enhances the effects. The NLChQM results behave like a mixture of the PV and 
PS (NJL) schemes. It can be understood by expanding quark-PS meson coupling term ~ Mk by its momentum 
k. 

• The fragmentation functions from NLChQM are relatively symmetric with respect to z whereas the local in- 
teraction models give asymmetric or tilted curves. After high-Q^ evolution, we find substantial deviations in 
comparison to the empirical data, in particular, in the small z region. The inclusion of the quark-jet contribution 
is expected to help to resolve this discrepancy. 

• The minus-type (valance) quark distribution function reproduces the empirical data qualitatively well for various 

values, while the plus-type one indicates sizable deviations at small x region. 

As a conclusion, the present nonlocal-interaction model, NLChQM provides very distinctive features, which have not 
been observed in usual models with local interactions, and contains interesting physical implications. We are working 
on the more realistic ingredients, such as the resonance and the quark-jet contributions. The result will appear 
elsewhere. The extension of the present theoretical framework to the flavor SU(3) octet mesons for the kaons is also 
under progress. 
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FIG. 5: (Color online) Minus- type distribution functions — /J) (left) and plus-type x(/^ + fu) (right) at different 

values, i.e. = 0.36 GeV^ and 4GeV^. The evolution is performed at LO and NLO. The empirical data are taken from 
Ref. g^. 




FIG. 6: (Color online) Valance-quark distribution function, x(/^ — fu) via the DGLAP evolution at LO to = 27GeV^ The 
empirical data are taken from Ref. [16j. 
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